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ABSTRACT 

This paper presents a search for radio transients at a frequency of 73.8 MHz (4 m wavelength) 
using the all-sky imaging capabilities of the Long Wavelength Demonstrator Array (LWDA). The 
LWDA was a 16-dipole phased array telescope, located on the site of the Very Large Array in 
New Mexico. The field of view of the individual dipoles was essentially the entire sky, and the 
number of dipoles was sufficiently small that a simple software correlator could be used to make 
all-sky images. From 2006 October to 2007 February, we conducted an all-sky transient search 
program, acquiring a total of 106 hr of data; the time sampling varied, being 5 minutes at the 
start of the program and improving to 2 minutes by the end of the program. We were able to 
detect solar flares, and in a special-purpose mode, radio reflections from ionized meteor trails 
during the 2006 Leonid meteor shower. We detected no transients originating outside of the solar 
system above a flux density limit of 500 Jy, equivalent to a limit of no more than about 10 -2 
events yr _1 deg -2 , having a pulse energy density > 1.5 x 10~ 20 J m 2 Hz -1 at 73.8 MHz for 
pulse widths of about 300 s. This event rate is comparable to that determined from previous 
all-sky transient searches, but at a lower frequency than most previous all-sky searches. We 
believe that the LWDA illustrates how an all-sky imaging mode could be a useful operational 
model for low-frequency instruments such as the Low Frequency Array, the Long Wavelength 
Array station, the low-frequency component of the Square Kilometre Array, and potentially the 
Lunar Radio Array. 

Subject headings: instrumentation: interferometers — methods: observational — radio continuum: gen- 
eral 

1. Introduction 
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Transient emissions — in the form of bursts, 
flares, and pulses from compact sources — are the 
signposts for explosive or dynamic events. As 
such, transient sources offer insight into a variety 
of fundamental aspects of physics and astronomy, 
ranging from studying the mechanisms of particle 
acceleration on the Sun and nearby stars to track- 
ing stellar evolution and death across the Universe 
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to probing the intervening medium(a) . 

At radio wavelengths, there are well-known 
classes of transients, such as the Sun and ra- 
dio pulsars, as well as a long history of observ- 
ing transients from triggers at other wavelengths, 
s uch as 7-ray burst (GRB) afterglows (e.g. , 



Pool ev et al.l 119971; iTavlor et al . 1998; Frail et al 



2000), or the monitoring of known transient 



sourc es such as X-ray binaries and microquasars 
e.g., IWaltman fe Ghigol [19981 iKlein-Wolt et al ' 



2002 iMartf et alj|200 j iRupen et alj|2002jl . Fur 

ther, a series of observations and discoveries 
over the past decade have emphasized that 
the radio sky may be quite dynamic. Known 
sources have been discovered to behave in new 
ways and what may be entirely new classes 
of sources have been discovered — pulsed ra- 
dio em ission has been ob served from brown 
dwarfs ( Hallinan et al. 2007t) and form erly radio- 
quiet magnetars ( Camilo et al.l 120061) ; single or 
highly intermittent pulses have been detected 
from neutron stars, als o known as rotating ra- 



dio transients (RRATs, McLaughlin et al. 2006 



iKeane et al1l2010t iBurke-Spolaor fc Bailed l2010l ) 
intense giant pulses have been det ected from the 
Crab pulsar (jHankins et al.l 120031 ) 
as-yet unidentified radio 



and several 
transients have been 



Lorimer et al 



20071) . 



found (iHvman et al.1 120051 iBower et al.1 12007 



There are also numerous classes of objects that, 
by extension of known physics, have been hypoth- 
esized to be radio emitting, and which might ap- 
pear as radio transients. T hese include extrasolar 
planets (|Farrell et al.lll999h . magnetar flares or gi- 
ant pulses from pulsars in other galaxies (e.g., 
McLaughlin fe Corded 120031). prompt emission 



from GRBs (|Usov fc KatzfeOOOtlSagiv fc Waxman 
2002 ). evaporating black hol es (iRees 1977), an d 



extraterrestrial transmitters ( Cordes et al.lll997 ). 

Meter-wavelength observations have a long his- 
tory of being used for studying the time domain. 
For example, solar radio emissions have long been 
known to be dynamic, and pulsars were first dis- 
cover ed in 81 MHz observations (jHewish et al 
1968). Advantages of low-frequency instruments 



for surveying the transient sky include naturally 
large instantaneous fields of view and sensitiv- 
ity to steep spectral index sources such as might 
result from coherent emission processes, while dis- 
advantages include radio- wave propagation effects 



such as dispersion and multi-path propagation, 
absorption either within the source or along the 
line of sight, ionospheric disturbances, and ra- 
dio frequency interference (RFI). A number of 
new telescopes with meter-wavelength capabili- 
ties are under construction, includ ing the Murchi- 
son Wide-field Array (MWA, lLonsdale et al.l 
20091 ). the Low Frequency Array (LOFAR), the 
Precision Array to Probe the Epoch of Reion- 
ization (PAPER, IParsons et"aH |2010|). and the 



Long Wavelength Array (LWA, lEllingson et al 
2009) — many of the science cases for these tele- 



scopes explicitly include searches for radio wave- 
length transients. It is anticipated that the suc- 
cess of these telescopes will motivate the low- 
frequency c omponent of the Squa re Kilometre Ar- 
ray (SKA, iDewdnev et al.ll200E) and the Lunar 



Radio Array (LRA, lLazio et al 



2009). 



This paper presents an all-sky monitoring cam- 
paign for radio wavelength transients conducted 
on the Long Wavelength Demonstrator Array 
(LWDA) . The objective is two- fold, both to detect 
or constrain radio transients from the observations 
and to illustrate how future instruments could be 
used for all-sky searching and monitoring obser- 
vations. The plan of this paper is as follows. In 
<j2] we describe the LWDA itself and its data ac- 
quisition path, in fj3]we describe the LWDA data 
acquisition and imaging pipeline specific to the 
formation of essentially all-sky images and radio 
transient searching, in §3] we present our results 
and analysis, in Jj5] we discuss possibilities for fu- 
ture instruments, and in f}6] we summarize our 
conclusions. 

2. The Long Wavelength Demonstrator 
Array 

The Long Wavelength Demonstrator Array 
(LWDA) was a testing platform developed by the 
Naval Research Laboratory prior to construction 
of the LWA. The LWDA site work included field- 
ing prototype hardware for the LWA, development 
of site preparation techniques, continual monitor- 
ing of the RFI environment, software develop- 
ment, and initial science. Located on the Plains 
of San Agustin in New Mexico near the center 
of the National Radio Astronomy Observatory's 
Very Large Array (VLA), the LWDA consisted 
of 16 dual-polarization dipole stands operating 



2 




as a phased di pole array with a frequency range 
of 60-80 MHz (I York et alJl2007i Figured ). The 
element locations represented a compact subset 
of a larger 256 stand pseudo-ran dom distribution 
that was developed for the LWA (|Kogan fe Cohen 
2009h . The maximum baseline for the LWDA 
was approximately 20 m (Figure [2]), providing 
an angular resolution at the zenith of about 12° 
at 74 MHz. 

Each of the 16 LWDA dipole stands consisted 
of 2 "droopy," fat dipoles, one per polarization, 
with the dipoles mounted at a 45° droop angle to 
a central post (Figure [3]) . Each dipole consisted 
of two "blades," each 28 cm wide by 1.05 m long. 
This design was chosen to increase the symme- 
try in the E- and f/-planes ( Kerkhoff fc Ellingsonl 
2005), broaden the antenna beam pattern (i.e., in- 
crease antenna sky coverage), and obtain a larger 
operational bandwidth. 

The radio frequency (RF) signals from each 



dipole were amplified by an active balun (24 dB 
of fixed gain) before being passed to the LWDA 
electronics, located within a central, shielded shel- 
ter. Inside the shelter, the incoming RF data 
stream entered the digital receiver signal process- 
ing chain. The signals were digitized by a dual- 
channel 10-bit analog-to-digital converter (ADC) 
at 100 megasamples per second. 

The LWDA electronics had the capability to 
provide digital delay beamforming of two fully in- 
dependent beams of bandwidth 1.6 MHz each as 
well as the capability for all-sky monitoring. This 
capability was implemented through two identi- 
cal signal processing chains that had independent 
control parameters to allow for two frequencies, 
two spatial beams, or two polarizations. The first 
stage of digital processing was a first-in-first-out 
(FIFO) pipe that provided an integer sample de- 
lay from ns to 150 ns in steps of 10 ns. This 
delay provided partial compensation of geometric 
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Fig. 2. — Location of the antennas within the 
LWDA. 



path and system time delays at this stage. Fol- 
lowing the FIFO was a complex mixer where the 
in-phase (I) and quadrature (Q) components for 
each signal chain were generated using a CORDIC 
rotator. The signals then passed through a low- 
pass filtering stage using a Cascaded Integrator 
Comb (CIC) filter and were decimated by a fac- 
tor of 14. Next, the signals passed through an- 
other FIFO, which provided coarse sample delays 
in excess of 2 /xs in 140 ns increments. The final 
1.6 MHz bandpass for each data stream was set 
in a finite impulse response (FIR) filter that also 
reduced the data rate by another factor of 3. Fol- 
lowing this stage, the data were interleaved and 
passed to an adder board at the final data rate of 
approximately 4.8 MSPS. The adder board was 
used to either sum (beamform) or interleave (all- 
sky image) the incoming coherent signals. 

The LWDA frequency range was chosen, in 
large part, because early plans for the array in- 
cluded the goal of operating it as o ne of the ele- 
ment s of the VLA 74 MHz system ( Kassim et al.l 



2007), which in turn was centered on the (primary) 
frequency allocation for radio astronomy at 73.0- 
74.6 MHz in the U.S. The operational frequency 
range was much broader and includes frequency 
allocations for a number of different services, pri- 
marily TV broadcasting, at the time of the LWDA 
operation!]] Tests with the LWDA and the pro- 




Fig. 3. — Photo of several LWDA elements show- 
ing the droopy blades and an NRAO VLA antenna 
in the background. 



totype LWA equipment have shown the presence 



this band (iDuffin & Ravi 2007; 


Hicks et al. 


2007; 


Paravastu et al. 2007; Rav et al. 


120071 Craid 


2008; 


Jaeger & Mute! not dated). Further, experience 



with the VLA 74 MHz system indicated that, for 
that system, the most significant source of inter- 
ference was the VLA itself. Indeed, in comparison 
with the VLA, for which the 74 MHz feeds were 
mounted near prime focus (several meters above 
the ground) , the LWDA should be less susceptible 
to RFI. Our experience with the LWDA was con- 
sistent with the experience from the VLA, in that 
no significant amount of observing time was lost 
due to RFI. (See also jgH ) 

The LWDA system was sky-noise dominated by 
at least 6 dB over its full operating range. The 
LWDA construction was completed and first-light 
images were obtained on 2006 October 23. 

Our focus here is on the all-sky observations 
obtained by the LWDA, but it also had a two- 
element interferometer mode. In this mode, an 
outlier dipole could be utilized to monitor the flux 
densities of isolated point sources, contributing 
to a st udy of the secular flux dens ity decrease of 
Cas A ()Helmboldt fc Kassim|[2009h . 



from analog to digital TV broadcasts, and much of the 
LWDA operating frequency range is no longer used for TV 
broadcasts at the time of writing (Crane 2009). 



1 Since the time of these experiments, the U.S. has converted 
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3. All-Sky Observations 



FITS-IDI format (Clar ke et al. 



2007). 



The observations presented here consist of op- 
erating the LWDA in its "all-sky" imaging mode 
observing at zenith. In this operational mode, 
the all-sky data acquisition system cycles serially 
through all 120 pairs of dipoles in the array and 
correlates the data from each pair. Each base- 
line is measured for 51 ms with a total cycle time 
of 13 s for both polarizations of all baselines (in- 
cluding auto-correlations). As a result of various 
processing overheads, the on-sky cadence of imag- 
ing described herein is slower. Our observations 
began with a 5 minute cadence (i.e., 13 s to acquire 
the full set of cross correlations to form an image 
with images repeated every 5 minutes), which was 
later improved to a 2 minute cadence for most of 
the data reported here. 

Our focus on an (all-sky) imaging transient 
search stems from the fact that the number of 
dipoles in the LWDA was relatively small so that 
cross-correlation of all of them was computation- 
ally tractable. The alternative to an imaging 
search, in which transients are identified by com- 
paring images of the sky at different epochs, is a 
non-imaging search, in which transients are iden- 
tified in time series of the voltages or intensities 
from a telescope (s). Traditionally, the tradeoffs 
between the two kinds of searches have been that 
imaging searches obtain higher angular resolution 
at the cost of time resolution with the converse 
being true for non-imaging searches. The choice 
between the kind of transient search to employ 
depends not only on the instrument available, but 
also upon what is kno wn about the t ransient pop- 
ulation^) of interest ( Nemiroff 20031) . 

The observations were centered at a frequency 
of 73.8 MHz with 80 spectral sub-bands over a to- 
tal bandwidth of 1.6 MHz (20 kHz per sub-band). 
During the course of the array installation, the 
amplitude and phase of the signal path from each 
dipole were determined. These values were applied 
in real time as the cross-correlations (visibilities) 
were formed. The normalized bandpass shapes 
were applied to correct for individual telescope 
variations. The visibilities were then recorded on 
disk for later processing. 

Post-processing consisted of the following steps: 



2. The visibility files were ingested by the 
Astronomical Image Processing System 
(AIPS)| 

3. Standard coplanar imaging procedures were 
used to convert the visibilities to images. 

4. Strong sources in the images were blanked 
(see below), and, because the images for- 
mally extend beyond the horizon, they were 
clipped at the horizon. 



Figure |4] presents a typical all-sky image. Sources 
easily detected in the LWDA all-sky maps include 
Cas A, Cyg A, the inner Galactic plane, North 
Polar Spur, Loop III, and the Sun. 

We used the apparent brightness of both Cas A 
and Cyg A to obtain a crude estimate of the power 
pattern of the LWDA. Simulati ons of the power 
pattern of the individual dipoles ( York et al.ll20071 ) 
suggest that the gain decreases toward the horizon 
(as would be expected from a dipole) as well as to- 
ward the zenith. The latter is caused by the design 
goal of operating the dipoles over a relatively large 
frequency range. The brightness measurements of 
Cas A and Cyg A do show these large-scale gain 
variations as well as smaller variations, which we 
estimate to be approximately a factor of 2 (3 dB) 
in amplitude. 

The all-sky transient data were passed into an 
analysis pipeline written within the AIPS software 
system. The pipeline first clipped high phase- 
center amplitudes that are known to be generated 
by sporodic monitor and control (M&C) software 
issues. The data were then Fourier transformed 
into dirty images of the entire sky visible at the 
LWDA site. Even though the LWDA can ob- 
serve essentially the entire sky, "wide-fiel d" imag- 
ing techniques (jCornwell fc PerlevHl992l ) are not 
required given the size of the LWDA, its opera- 
tional frequency, and our use of snapshot imaging. 

Images were blanked around the positions of 
the strong sources that could be detected easily 
(i.e., Cyg A, Cas A, the Sun, the inner Galac- 
tic plane, Loop III, and the North Polar Spur). 
This strategy of blanking the images allowed us 
to use a relatively simple threshold test to search 



1. The raw visibilities were converted into the 
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Compariso n of all-sky LWDA images (green contours) with the 408 MHz all-sky image (color) from 



lHaslam et al.l (|l982t) . Insets in both panels (demarcated by the white circle) show the all-sky LWDA images 
(in color); in these insets, north to up and east is to the left. The dotted white line indicates the horizon for 
the LWDA, with the portion above the dotted white line being visible. The LWDA images are at 73.8 MHz 
and were acquired on 2006 October 28. Top: a snapshot acquired when the Galactic plane was passing 
nearly overhead. Individual, well known sources and the inner Galactic plane dominate the image. Bottom: 
a snapshot acquired when the inner Galactic plane was below the horizon. While the LWDA image would 
appear to be only noise, there are features that can be identified with the North Polar Spur and Loop III 
radio structures. 
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for the presence of transients, at the cost of reduc- 
ing somewhat our sky coverage. Near the end of 
the LWDA's operation, we did attempt to change 
the array so that the time sampling was in sidereal 
time. Had this been successful, it would then have 
been possible to subtract images acquired on dif- 
ferent days and construct "difference" images for 
which this blanking scheme probably would not 
have been necessary. 

As a measure of our ability to detect a tran- 
sient, particularly in light of our non-simultaneous 
sampling of the baselines, we conducted a num- 
ber of simulations in which transients of differ- 
ent strengths and positions within the image were 
inserted into the data. Table []] summarizes the 
results: Column 1 indicates what fraction of the 
baselines are illuminated by the simulated tran- 
sient, Column 2 shows the implied duration of 
the transient, and Column 3 shows the approx- 
imate signal-to-noise ratio (S/N) that the tran- 
sient would have to have in order to be detected. 
One important aspect of these simulations is that 
the simulated transients were inserted into the vis- 
ibility data. As such, these simulations do not 
take into account the antenna power pattern. The 
trend in Table Q] of increasing S /N required in or- 
der to detect the transient is perhaps not surpris- 
ing, but we find it reassuring that we could detect 
transients, even those whose durations are only a 
minor fraction of the time it took the LWDA to 
acquire a full sample of visibility data. 

Using our simulations and visual inspection of 
the images as a guide, each all-sky snapshot image 
was searched for events above the 5a level. Any 
image containing a potential transient event was 
saved for further analysis. 

4. Results 

We begin this section by both motivating and 
illustrating the capability of the LWDA to detect 
radio transients, and we then turn to the detection 
of astronomical radio transients. 

4.1. Detected Radio Transients: The Sun 
and Meteor Trails 

The Sun is a well-k nown variable radio sourc e 
at tens of MHz (e.g., IWild fc McCreadvl Il950ft . 
While the Sun is currently in a state of excep- 
tionally low radio emission, during the course of 



the LWDA observations reported here, a number 
of solar radio bursts were detected, including a se- 
ries of intense radio bursts on 2006 December 14 
(Figure [5]), probably associated with the X-class 
flares that occurred then. 

The focus of our monitoring program was pri- 
marily on possible extrasolar system transients, 
and, as such, we made no effort to track the solar 
flux density (32]). However, a similar monitoring 
program could be used to track the flux density 
of the Sun at these frequencies. Such a program 
would be more valuable if the array monitoring 
the Sun had both a larger bandwidth and a higher 
time resolution. With these improved characteris- 
tics, not only would such an array be capable of 
detecting solar radio bursts, it might be able to 
image their evolution in frequency. 

Conversely, as Figure [5] illustrates, the Sun can 
become a source of "interference." During these 
intense solar radio bursts, the dynamic range of 
the LWDA and imaging pipeline were sufficiently 
limited that only the Sun was visible. Whether 
other observations can be conducted during so- 
lar radio bursts will depend upon the imaging dy- 
namic range of future systems, however, exceeding 
an imaging dynamic range of order 100 seems dif- 
ficult for rapid, all-sky imaging pipelines. 

Another well-known class of radio transients 
is reflections from ionized meteor trails. As a 
meteor enters Earth's atmosphere, its velocity is 
high enough that it can produce an ionized trail, 
which in turn can refl ect a radio transmission (e.g., 
iMillman et~aDll948h . 

In 2006 November, during the Leonid meteor 
shower, data were acquired with the LWDA tuned 
to a frequency of 61 MHz. At the time of these ob- 
servations, this frequency was within the band al- 
located for television broadcasting, and there were 
a number of TV stations in New Mexico and the 
southwest U.S. that used this allocation. It is a 
simple matter to show that even a TV transmit- 
ter of relatively modest power located more than 
1000 km distant would still produce a strong re- 
flection. For the purposes of this observation, the 
process for data acquisition was modified to re- 
duce the number of bits and allow simultaneous 
sampling of all 16 dipoles. Figure [6] illustrates 
two examples of a meteor reflection imaged by the 
LWDA. 
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Table 1 

Transient Detection Simulations 



Baseline Fraction 


Time 


S/N 


Illuminated 


(s) 




100% 


6 


5 


90% 


5.4 


7 


75% 


4.5 


10 


50% 


3 


13 


30% 


1.8 


30 



Note. — The results of injecting 
simulated transients into the LWDA 
data. Column 1 indicates what frac- 
tion of the baselines are illuminated 
by the simulated transient. Because 
the LWDA acquired all-sky imag- 
ing data by sampling baselines se- 
quentially, Column 2 shows the im- 
plied duration of the transient. Col- 
umn 3 shows the approximate signal- 
to- noise ratio (S/N) that the transient 
would have to have in order to be de- 
tected. 




Fig. 5. — Illustration of the LWDA's detection of a solar radio burst. Left: the sky before the solar radio 
burst. Right: the sky during the solar radio burst of 2006 December 14. This example is one of the several 
intense radio bursts that occurred in 2006 December. 
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Fig. 6. — LWDA images of reflected TV signals from ionized meteor trails during the 2006 Leonid meteor 
shower. Each image is constructed from a 0.1 s data acquisition at a frequency of 61 MHz. Top left and top 
right: images of reflected signals from Leonid meteors. Bottom: a reference image showing the nominal sky 
observed by the LWDA. The zenith is at the center of the image. The two discrete emitters near the center 
of the field are celestial radio sources (Cas A and Cyg A). The band stretching from the upper right to lower 
left, terminating in the strong emission near the edge of the field is from the Milky Way. 
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4.2. Astronomical Radio Transients 

The observing campaign lasted from 2006 Oc- 
tober 27 to 2007 February 17, during which there 
were 59 days on which some useful data were ac- 
quired. On 2006 December 13, we were able to 
improve the cadence at which images were formed 
from 5 minutes to 2 minutes. The resulting to- 
tal observing time for the transient campaign was 
106 hr. During these observations, no transients 
were detected above the 5a level in any image. 
Based on comparative analysis of the brightnesses 
of Cas A and the Sun, we estimate the (la) noise 
level in the LWDA all-sky snapshots to be roughly 
500 Jy. We also assume that the effective sky 
coverage of the LWDA is roughly 10,000 deg 2 
(~ 7r sr). Although, we show essentially images 
of the entire hemisphere, the gain of the LWDA 
was unlikely to be constant across the sky, an as- 
sumption partially supported by modeling. Conse- 
quently, we assume that the effective sky coverage 
is somewhat less than a full hemisphere. 

A total of 29,437 data sets were processed 
through the imaging pipeline. Of these, 1764 pro- 
duced potential transient candidates. Visual in- 
spection of these showed that they could be ex- 
plained as non-astronomical transients: slight er- 
rors in our calculations of the extents of strong 
sources resulted in one pixel of a strong source 
not being blanked appropriately. 

Combining these results, we place a limit of 
< 10 -2 events yr _1 dcg~ 2 , having a pulse energy 
density > 1.5 x 10~ 20 J m" 2 Hz" 1 at 73.8 MHz 
for pulse widths of about 300 s. We now con- 
sider multiple characterizations of this result for 
easy comparison with existing results in the liter- 
ature. Our focus is on compari son with MOTOR 
(|Amv et al.ll 19891 ) and STARE (jKatz et al.ll2003l ). 
as both represent nearly all-sky surveys at frequen- 
cies below 1 GHz. 



1 pc 



10 kpc 



Following iKatz et all (|2003l) . the LWDA flux 
density limit can be recast in terms of the bright- 
ness temperature that an object would have to 
have in order to be detectable. Using parameters 
appropriate for the LWDA, and assuming a quite 
conservative 5000 Jy flux density limit (10er), we 
find 



T B > 40 K 



(1) 




20 

Log[Distance (cm)] 



for an object of linear size L located at a dis- 



Fig. 7. — Brightness temperature limits for ob- 
jects of varying sizes (diagonal lines) as a function 
of the distance to the object. Lines for fiducial 
object dimensions are shown as are fiducial dis- 
tances. The required brightness temperature for 
an object to be detectable is above and to the left 
of a line of constant object dimension. The range 
of the ordinate is approximately that allowed for 
the LWDA observations (viz. Figure [5]). 

tance D. Figure [7] shows the resulting bright- 
ness temperature limits for transients detectable 
by the LWDA. Sufficiently strong bursts from a 
neutron star or a "super-flare" from a nearby star 
are among the possible sources that might have 
been detectable by the LWDA (or an instrument 
like the LWDA). We consider these possibilities 
further in fj5] 

Table [2] compares our rate with those of other 
"all-sky" surveys at other wavelengths from the 
literature. In constructing Table [2 we have 
restricted our focus to other, largely similar 
search programs, namely those at frequencies be- 
low 1 GHz using dipole-based arrays. We have 
also restricted Table[2]to programs that conducted 
blind surveys, as opposed to those that targeted 
known sources (e.g., pulsars, GRBs). 

In general, our LWDA search program produces 
an upper limit to the event rate that is compara- 
ble to or competitive with other search programs, 
typically at higher frequencies. However, our se- 
rial data acquisition scheme results in a somewhat 
poorer limit on the pulse energy density. If the 
LWDA data acquisition system had had a more 
rapid sampling time, e.g., on a scale of order 0.1 s 
(the time scale for the acquisition of a signal from 
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Table 2 
Transient Event Rates 



Pulse Scaled 



Program 


Rate Limit 
(events yr _1 deg -2 ) 


Energy Density 
(J m" 2 Hz" 1 ) 


Timescale 
(s) 


Frequency 
(MHz) 


Pulse Energy Density 
(J m~ 2 Hz^ 1 ) 


Reference 


LWDA 


10~ 2 


1.5 x 10" 20 


300 


73.8 




this work 


HM74 


0.05 


~ 10" 23 


~ 0.1 


270 


3.4 x 10~ 23 


1 


K+77 


4 x 10~ 3 


2.5 x 10" 22 


> 1 


370-550 


1.3 x 10~ 21 


2 




4 x 10~ 3 


3 x 10~ 22 


> 1 


38-60 


1.5 x 10" 22 


2 


MOTOR 


163 


1Q -28 


0.025 


843 


1.1 x 10" 27 


3 


STARE 


10" 3 


1Q -24 


0.125 


611 


8.3 x 10" 24 


4 



Note. — The scaled pulse energy density is t he pulse energy den s ity sc aled to a frequency of 73.8 MHz assuming a 
nomin al sp ectral index of — 1. References: (1) iHuguenin fc Moore] |l974|); (2) iKardashev et al.l (|l977l ); (3) lAmv et ail 
(|l989l ): CD lKatz et aT] (|2003h . 



an individual dipole), our pulse energy density 
limit would have been comparable to many of the 
other search programs. 

An alternate means of characterizing radio 
transient searches is in terms of a "phase space" 
diagram (Co rdes et al. 2004b), which can be cast 
in terms of the brightness temperature of a radio 
transient in the Rayleigh- Jeans limit (Equation[TJ . 
Figure[5]shows the region of phase space probed by 
the LWDA in comparison with the range defined 
by various known and hypothesized phenomena. 
In producing this plot, we have assumed a lower 
distance limit to any celestial radio transient to 
be 1 pc, such as might be expected from a nearby 
low-mass star, brown dwarf, or planet ([JSJ. In 
practice, the local density of sources may be suffi- 
ciently low that there are no transients within 1 pc, 
which would have the effect of increasing the min- 
imum value of SD 2 that the LWDA could have 
probed. Clearly the upper distance limit could be 
the edge of the observable Universe, which would 
imply correspondingly higher brightness tempera- 
tures. 

As a final approach to constraining the event 
rate for transients, we consider a probabilistic one. 
The typical determination of an event rate (such 
as those in Table [2]) is derived from the total ob- 
serving duration T and the solid angle il observed. 
Instead, we ask what is the probability p of detect- 
ing a transient in a given observation, given that 



no transients were detected in our series of obser- 
vations. 

We conducted 29,437 trial observations dur- 
ing the c o urse o f this observational program. As 
DeGroot ( 19861 ) discusses, the probability of ob- 
taining i detections in N trials, given that the de- 
tection probability in any individual trial is p, can 
be related to the beta distribution 



P(n,j) 



T{n + 2) 



r(j + i)r(n-j + i) 



(2) 



where T(x) is the gamma function, j = 53i=i Xi, 
and Xi = 1 for a detection or for a non-detection. 
That is, j is simply the count of the total number 
of detections. 

With j = 0, and integrating Equation ([2]) to 
form the incomplete beta function, we can deter- 
mine the maximum detection probability p for a 
specified confidence level. For the LWDA or an 
array with LWDA-like parameters (« 10,000 deg 2 
sky coverage, 73.8 MHz frequency, < 300 s du- 
ration), we find that the probability of detecting 
a transient in a single LWDA observation cannot 
exceed p = 2 x 10~ 4 lest the probability that we 
should have detected a transient in the full LWDA 
transient campaign exceeds 99.7%. 

Given the relatively large solid angle accessed 
by the LWDA, only a modest improvement in 
the probability of detecting transients would re- 
sult from increasing the sky coverage for future in- 
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-9-8-7-6-5-4-3-2-1 1 23456789 10 11 12 



log l/W (GHz s) 

Fig. 8. — Phase space defined by radio transients. The abscissa is an uncertainty-like quantity given by the 
product of the observing frequency v and the duration of the transient W. For our LWDA survey, we define 
two regions. The first is from 2 to 6 s and is the interval required to accumulate the visibility data from all of 
the interferometric baselines ([J3]and Table [lj. The second is from 5 minutes to 24 hr and covers the cadence 
at which individual images were made to the longest continuous observing duration during our program. 
The ordinate is a pseudo- luminosity S p kD 2 , for a peak flux density S p k and distance D. We take a lower 
distance limit of 1 pc, such as might be expected from a nearby low-mass star, brown dwarf, or extrasolar 
planet. Various classes of known or hypothesized transients are indicated. The sloping solid lines indicate 
lines of constant brightness temperature, with a brightness temperature of 10 12 K taken as the dividing line 
between coherent and incoherent emission. 
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struments (or would have resulted from increasing 
the sky coverage of the LWDA itself). Rather, an 
increased probability is likely to result only from 
improved sensitivity, developing a search method 
that does not require blanking the Galactic plane, 
or both. 

5. Future Possibilities 

In the near future, the first station of the Long 
Wavelength Array (LWA-1), the LOFAR, and the 
International LOFAR Telescope (ILT) will become 
operational, with LOFAR being a subset of the 
ILT. All will consist of dipole stations having an 
operational frequency range overlapping or com- 
parable to that of the LWDA. 

The individual dipole stations of LOFAR and 
the ILT will have a factor of a few more dipoles 
than the LWDA while LWA-1 will 256 dipoles 
(16 x as many), but they may not be able to im- 
prove upon our flux density limits substantially 
because the diameters of the LWDA, LWA-1, and 
LOFAR stations are sufficiently small that they 
are confusion limited on short time scales. How- 
ever, with multiple LOFAR stations, the identifi- 
cation and excision of RFI m ay be improved sub- 
stantially (jBhat et al.l 120051) . In the case of the 
LWA-1, all-sky imaging will be carried out by pass- 
ing a narrow bandwidth (~ 100 kHz) signal from 
each of the 256 stands to a dedicated cluster for 
correlation and imaging. This observing mode op- 
erates in tandem with the beamforming modes so 
that it will be possible to image the visible hemi- 
sphere full time. 

The MWA will also conduct transient searches, 
at a somewhat higher frequency range than the 
LWDA, but it will likely not include an all-sky 
imaging capability. Building upon these next- 
generation low-frequency arrays will be the SKA 
and the LRA. Their designs will be influenced by 
the work on LOFAR, MWA, and similar low ra- 
dio frequency interferometers, but, in both cases, 
we anticipate that transient searches, and poten- 
tially all-sky imaging, will be a capability. As Fig- 
ure [5] illustrates, and as we now discuss below, in 
principle, an imaging dipole station, similar to the 
LWDA, could probe a variety of known or hypoth- 
esized classes of transients. 



5.1. Radio Pulses from Ultra-high Energy 
Cosmic Rays or Neutrinos 

Intense, short-duration pulses (~ 1 MJy in 
~ 10 ns) at decimeter and meter wavelengths 
have been detected from the impact of ultra-high 
energy cosmic rays on Earth's atmosphere (e.g., 
^' l965llWeekell200lL IFalcke fc Gorhaml 
least one case even imaged 
not dissimilar to the LWDA 
20051) . High-energy neutrinos im- 



Jelley et al 
20031 ). 



by 



and 
radio 



m at 
a rray 



( Falcke et al 



pacting the lunar regolith should also produce ra- 

dio p u lses via the Askar'yan effe ct ( Dagkesamanskii fe Zheleznvkhl 



1991 



19891 : IFalcke fc Gorhaml 12003). though no such 
pulses have been d etected_| from the Moon to 
date (jHankins et al 
Bercsnv ak et al.l2 005; 



Gorham et all 12004 



Langston et al. 2009; Scholtcn et al 



2009tlBuitink et al.l2010l : ljames et al-hoidljaeger et al 
120101 ). In the case of the lunar neutrino ra- 
dio pulses, the expected pulse amplitude is more 
poorly constrained because it depends, in part, 
upon properties (pri marily surface roughness) 
of the lunar surface jAlvarez-Muniz et al. 
Scholten et al. 20061 : Ijames fc ProtheroeT 



Nonetheless, if sufficiently energetic neutrinos ex- 
ist, it is plausible that radio pulses would be emit- 
ted at frequencies relevant for an LWDA-like in- 
strument. 

The intrinsic width of these radio pulses is ex- 
tremely short, far shorter than the 51 ms data 
acquisition time for an individual baseline of the 
LWDA; as an example, IFalcke et al.l (|2005l ) mea- 
sure a width narrower than 30 ns. As noted above 
(<j3]), during the commissioning of the LWDA, a 
sporadic M&C issue was recognized in that occa- 
sionally a baseline would produce an anomalously 
large amplitude. Upon inspection, we found that 
baselines with anomalously large amplitudes also 
had a phase of 0° (which is one of the reasons 
that we explain these anomalous amplitudes as 
being due to an M&C issue). Furthe r, unlike the 



LOFAR Prototype Station (LOPES, IFalcke et al 



2005), we had no co-located particle detector to 
use as a trigger for cosmic ray detection. Subse- 
quently, we employed a simple threshold test, in 
which baselines with anomalous amplitudes and 
phases near zero were excised. A consequence of 



3 No pulses from the Moon have been detected, but the 
Askar'yan effect has been demonstrated in terrestrial ac- 



celerators l lGorham et al.| [2005. 2007). 
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this flagging methodology is that radio pulses from 
high-energy particles at the zenith would not be 
recognized by our analysis of LWDA data. 

5.2. Stars and Substellar Objects 

The Sun generates intense emission at LWDA 
frequencies, notably Type III and IV radio bursts, 
but even the strongest such radio bursts are far too 
faint to be detected over interstellar distances with 
the current most sensi tive meter- or decameter- 
wavelength telescopes. Gergelv ( 1986h considered 
the detection of solar-type stars at low radio fre- 
quencies and finds that millijansky sensitivities 
would likely be required to detect the equiva- 
lent of the most intense solar radio bursts from 
nearby stars. There are approximately 20 solar- 
type stars in the sol ar neighborhood (< 10 pc, 
Holmberg et~aD l2009h . The absence of any tran- 
sients that could be identified with these stars puts 
only modest constraints on their levels of activ- 
ity. Not only is the LWDA sensitivity likely to 
be insufficient, if these stars have solar-like cycles 
over decadal time scales, the absence of any flar- 
ing could merely reflect that a star is currently in 
a quiescent phase of the cycle (much like the Sun 
is in the current solar cycle). 

Late-type stars (M dwarfs) and brown dwarfs 
also exhibit flaring activity, particularly at centimeter- 



Giidcl 



and meter wavelengths (> 300 MHz, e.g. 
20021; lHallinan et alj|2007h . Little is known about 
the emission from brown dwarfs at LWDA fre- 
quencies. However, the most sensitive searches 
for decameter-wavelength emission from late-type 
stars have, at best, marginal detections, even when 
those stars were observed t o be fla ring at shorter 
wavelengths iJackson et al. Even these 



marginally detected flares have flux densities of 
order 1 Jy. While more numerous and likely to be 
closer than solar-type stars, the LWDA sensitivity 
also means that we can place only modest limits 
on the likelihood of strong flares from nearby M 
dwarfs. 

Finally, within the solar system, in addition to 
solar radio bursts, Jupiter also produces intense 
radio bursts due to an ele ctron cyclotron maser 
within its magnetosphere ([Melrose! I2005T ) . This 
emission cuts off above around 40 MHz, a fre- 
quency determined by the ratio between the cy- 
clotron and plasma frequencies as a function of 
altitude within the Jovian magnetosphere. 



The magnetospheric emissions from solar sys- 
tem planets and the discovery of extrasolar plan- 
ets have motivated a number of both theoretical 



dZarka et al.lll997t iFarrell et al.lll999HZarka et al 
| 200l | :lFarrell et alj|2004 lLaziq et alj2004IStevens 



2005: Gricssmeier et al 



Griessmeier et al. 2007b 



2005] 
"3> 



work dYantis et all 119771 IWindee et all Il986t 



Zarkal 120061 l2007t 

and observational 



Bastian et al.l2000ULazio et al.ll2004l:lRvabov et al 
2004 iGeorge fc Stevens! l2007t lL azio fe FarrelJ 



2007t ISmith et al.U "2009) on magnetospheric emis 



sions from extrasolar planets, including some be- 
fore the confirmed discovery of any extrasolar 
planets. Unfortunately, even the most optimistic 
predictions for extrasolar planetary radio emission 
do not predict flux densities in excess of 1 Jy at fre- 
quencies near 100 MHz. Moreover, the most sen- 
sitive searches for extrasolar planetary radio emis- 
sion near 74 M Hz place sub-Jansky c o nstraints on 
this emission ( Lazio fc Farrefll l2007t E azio et al 
20091) . 



5.3. Pulsar Giant Pulses 

We consider giant pulses from radio pulsars 
as an exemplar of a short-duration, intense ra- 
dio transient that, in principle, could be de- 
tected throughout th e Local Group by the LWDA 
(jCordes et al.ll2004al ). For the specific case of the 
Crab pulsar, giant pulses h ave been seen at fre- 
quencies as low as 74 MHz dComella et al. 196£ ; 



Rankin et ail Il970h . and ISallmen et al] fl999) 



showed that the pulses are indeed broadband, at 
least over the frequency range 610-4900 MHz. 

The spectral index for Crab pulsar giant pulses 
is steep (a ~ —3.4, S v oc v a ), and current ob- 
servations fin d no upper limit to t he amplitude of 
giant pulses ( Sallmen et al. 1999t Lundgren et al 



19951 ) . Howe ver, by comparison to lo wer frequency 
observations. iLundgren et al.l (|l995l ) also find that 
the rate of giant pulses decreases with decreas- 
ing observational frequency, such that giant pulses 
at 800 MHz are 400 times more frequent than 
at 146 MHz. Scaling both the rate of giant pulses 
and the pulse amplitudes to the LWDA frequency 
of 73.8 MHz, we estimate that, on average, the 
Crab pulsar will produce a giant pulse once every 
50 hr with a flux density exceeding 1 MJy (i.e., 
- 2000cr for the LWDA, Tabled]). 

At higher frequencies, Crab pulsar giant pulses 
have extremely narrow widths, often being unre- 
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solved in time. With a pulse period of approxi- 
mately 33 ms, even if the Crab pulsar or a pul- 
sar like it produced a sufficiently strong pulse to 
be detected, only a single baseline would appear 
to be illuminated by our data acquisition proce- 
dure (51 ms per baseline, Sj3j) . However, at LWDA 
frequencies, pulse broadening becomes significant. 
Pulse broadening increases the duration of the 
pulse, at the cost of decreasing its amplitude. 

In addition to being frequency dependent, 
pulse broadening also depends upon direction 
through the Galaxy. Using the NE2001 model 
(jCordes fc Lazioll2002h . we estimate that the mag- 
nitude of pulse broadening can exceed 500 ms 
(~ 10 LWDA baselines being illuminated) for dis- 
tances of order 5 kpc through the Galactic disk, 
implying Galactic latitudes \b\ < 5°. Unfortu- 
nately, these low Galactic latitudes constitute only 
a small fraction of the total sky (< 1%), and we 
also blanked the Galactic plane as part of our ef- 
forts to remove strong sources from the images be- 
fore conducting the statistical tests. Thus, we con- 
clude that while giant-pulse emitting pulsars seem 
capable of producing sufficiently strong pulses to 
be detectable, either our serial data acquisition 
procedure would have not been sensitive to such 
pulses or significant pulse broadening would have 
rendered them too faint to be detectable. 

5.4. Gamma-ray Bursts (GRBs) 

We conclude this section by considering radio 
pulses from GRBs, as an exemplar of possible ra- 
dio pu lses of extragalactic ori g in. Usov fc Kat3 
(|200dh and ISagiv fc Waxmanl (|20o3T both have 
predicted that GRBs should have associated 
prompt emission, most likely below 100 MHz. Pre- 
dicted flux densities are highly model dependent 
and range from essentially undetectable to in ex- 
cess of 1 MJy. 

Observationally, there have been a number of 
searches for radio puls e s associated with G RBs. 
Cortiglioni et all (|l98lh . Ilnzani et al.l (jl982t ). and 
Amy et al.l (|1989) all detected some dispersed ra- 
dio pulses, but found no c onvincing associations 
with GRBs. iBalsanol (|l999h found a dispersed ra- 
dio pulse apparently coincident with GRB 980329, 
but it was narrowband, which has led to it be- 
ing interpreted as due to terrestrial interference. 
Various searches for radio pulses associated with 
GRBs (including precursor pulses) have been con- 



ducted at 151 MHz 



Koranvi et all (|l994 Il995l) 



Dessenne et al.l ( 1996 ). Typical upper limits have 



been approximately 100 Jy. 

5.5. Exotica 

One of the motivations for this, and similar, 
searches is simply that previously unknown classes 
of sources may be discovered. While it is difficult 
to assess the probability of an unknown class of 
sour ces, we can consider extrate rrestrial transmit- 
ters ( Cocconi fc Morrisonlll959l) as an example of 
an exotic population of sources 

As Loeb fc Zaldarriaga ( 20071) discuss, much 
of the human-generated radio radiation from the 
Earth is emitted in the range 50-400 MHz. In- 
deed, the upper end of the LWDA's operational 
frequency range was chosen to avoid FM ra- 
dio broadcasts, and 34.11 illustrates the detection 
of TV signals with the L WDA. In contrast to 
Loeb fc Zaldarriaga! (|2007l ). however, we consider 
short duration pulses, such as might originate from 
a transmitter on a planet rotating into or out of 
view, whereas they consider detecting signals in 
long integrations (~ 1 month). 

The LWDA's frequency sub-bands of 20 kHz are 
well matched to the typical bandwidths, at least 
for terrestrial transmitters in this frequency range. 
There is an increasing usage of so-called spread 
spectrum transmitters, with bandwidths (much) 
larger than 20 kHz. However, these are often at 
higher frequencies, and an LWDA-like instrument 
would still receive all of the power within its re- 
ceived bandwidth. 

Unfortunately, with a 500 Jy sensitivity, it is a 
relatively simple matter to show that a transmit- 
ter at a distance of 10 pc would have to have an 
effective isotropic radiated power (EIRP) of about 
10 17 W (170 dBW) to have been detectable with 
the LWDA. By contrast, the EIRP of even pow- 
erful military radars, such as the Air Force Space 
Surveillance System (AFSSS), only approach 5 x 
10 13 W (137 dBW). Further, given that all (or 
most) of the stars within 10 pc are cataloged, it 
is not clear that an all-sky search, such as the one 
pursued here, is the most effective search strat- 
egy, unless one postulates a set of freely floating 



Doing so is also appropriate given that this year is the 
50 th anniversary of the first search for radio transmissions. 
F. Drake's Project OZMA. 
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"beacons." 

5.6. General Considerations 

Examination of Figure [8] shows that an LWDA- 
like instrument might be able to probe a significant 
region in this phase space, and previous sections 
have discussed a range of source classes that the 
LWDA might have been able to detect, in princi- 
ple, even though the implementation of the data 
acquisition system resulted in no detections. In 
addition to the aforementioned issue that a more 
rapid time sampling would likely be required, are 
there other considerations for a dipole phased ar- 
ray? 

Much of the phase space accessible, or poten- 
tially accessible, to an LWDA-like instrument 
is for coherent transients, nominally taken to 
be those with brightness temperatures in excess 
of 10 12 K. Such transients will n ecessa rily evolve 
rapidly. ICordes fc McLaughlin ( 20031 ) have dis- 
cussed searches for fast radio transients. A key 
consideration for such searches is the effective time 
resolution, which is determined not only by the in- 
strumental characteristics (receiver bandpass) but 
also by dispersion smearing and radio-wave scat- 
tering. 

Following Cordes fc McLaughlin ( 20031 ). we 
consider two time scales related to the time reso- 
lution, and implied processing, for searches with 
an LWDA-like instrument. The highest time reso- 
lution that the instrument would be able to obtain 
is simply 

At Al/ ~ (Az/p 1 , (3) 

for a receiver frequency channelization or sub- 
band of Av, which for the LWDA (Av = 20 kHz) 
was 50 /is. Because of the dispersion relation of 
the interstellar medium, probing to any distance 
for fast pulses requires de-dispersion. As an illus- 
tration, we consider how far one could probe, ig- 
noring dispersion. We take as a criterion that the 
dispersion across the receiver sub-band produces a 
smearing comparable to ti me resolution implied by 
Equat ion ((3|). Then, from Cordes fc McLaughlin! 
(120031) . 



SDM = 120pccm~ 3 (A^ Hz 



-2, ,3 



^MHz; 



(4) 



of MHz, and the error in the dispersion mea- 
sure (5DM is in the canonical units of pc cm~ 3 . 
Alternately, <5DM can be considered the incre- 
ment in dispersion measure (DM) for searches 
that incorporate searching through dispersion. 
For the LWDA, SDM = 0.1 pc cm- 3 , imply- 
ing that a search for fast transients with the 
LWDA that did not incorporate de-dispersion 
could only probe to an effective DM of this value. 
For reference, PSR J0108-1431 has the lowest 
known DM, 2.38 pc cm" 3 , at a distance of 240 pc 



where Aj/kHz is the receiver sub-band in units of 
kHz, and ^mhz is the central frequency in units 



([Manchester et alJl20Q5T: iDeller et alj|2009h . 

Conducting a search that incorporated compen- 
sation for dispersion smearing would effectively 
mean conducting a de-dispersion search on each 
resolution elem ent (beam) within the field of view 
(ICordesll2009l) . In contrast to proposed searches 
at higher frequencies (~ 1 GHz), the much lower 
resolution at LWDA frequencies makes this a less 
onerous task. For instance, with the LWDA's res- 
olution of 12°, there are only about 150 indepen- 
dent resolution elements in the entire sky. Further, 
searching each is of course an embarrassingly par- 
allel problem. 

6. Conclusions 

We have described the Long Wavelength Demon- 
strator Array (LWDA) and its operation as an 
all-sky transient monitor. A 16-element dipole 
phased array, operating over the frequency range 
60-80 MHz, the LWDA was used as a technical 
test bed for the LWA. We have described the sig- 
nal flow from the antennas to the receivers and to 
a simple software correlator as well as the config- 
uration of the dipoles. 

The individual dipoles of the LWDA had a field 
of view comparable to the entire sky. In late 2006 
and early 2007, we used the LWDA as a tran- 
sient search instrument by making a series of all- 
sky images at frequencies of 61 and 73.8 MHz. 
The 61 MHz observations were designed explicitly 
to search for reflections from ionized meteor trail 
during the Leonid meteor shower in 2006. The 
73.8 MHz observations were designed to search for 
astronomical transients, and we acquired a total 
of 106 hr of data, with a time sampling ranging 
from 2 to 5 minutes between images. We were able 
to detect solar flares and, by utilizing a special- 
purpose mode, TV reflections off ionized meteor 
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trails during the 2006 Leonid meteor shower. We 
detected no radio transients outside of the so- 
lar system above a flux density limit of 500 Jy, 
equivalent to a limit of no more than about 10~ 2 
events yr _1 deg -2 , having a pulse energy density 
> 1.5 x 10~ 20 J in" 2 Hz" 1 at 73.8 MHz for pulse 
widths of about 300 s. This event rate is compa- 
rable to many existing limits from previous all-sky 
surveys, but at a lower frequency than most pre- 
vious all-sky searches. 

There are a number of emerging arrays (e.g., 
MWA, LOFAR, LWA) for which all-sky, or at least 
very wide-field, imaging will be a capability. All of 
these arrays will offer at least two advantages with 
respect to the LWDA. First, all are anticipated 
to be operational for a far longer time than was 
the LWDA, which should improve upon the event 
rates substantially. Second, all of these arrays will 
have a higher angular resolution imaging capabil- 
ity, so that if a transient were to be detected, much 
higher precision astrometric information could be 
obtained. 
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port in establishing the LWDA near the VLA. We 
thank K. Weiler, D. Munton, and L. J Rickard 
for their guidance during various stages of the 
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software to generate Figure HI This research has 
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Lunar Science Institute (via Cooperative Agree- 
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